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ABSTRACT 


Data  on  small  local  earthquakes  were  obtained  from  a  tripartite  located  on 
the  island  of  Honshu,  Japan.  Epicenters  and  origin  times  were  computed  from  the 
data  but  local  magnitudes  could  not  be  determined.  LASA  beams  were  formed  at  and 
around  the  computed  epicenters  to  determine  if  the  event  could  be  detected  on  the 
array.  From  the  original  population  of  about  150  events,  only  28  could  be  detected 
or  were  marginal.  The  LASA  epicenters  were  obtained  by  beamsplitting  and  were 
found  to  be  consistent  with  the  theoretical  location  errors  expected  for  beams  steered 
to  Japan.  Although  the  total  number  of  events  visible  at  LASA  were  small,  it  was 
consistent  with  a  beam  detection  threshold  of  3.  5  for  beams  steered  to  Japan. 
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I.  SUMMARY  OF  THE  EXPERIMENT 


At  the  time  of  the  installation  and  initial  operation  of  the  Large  Aperture  Seismic 
Array  (LASA)  in  Montana,  evaluation  of  its  performance  was  hampered  by  the  lack  of 
a  reference  population  of  small  teleseismic  events.  In  an  early  attempt  to  provide  such 
a  reference,  for  at  least  one  small  epicentral  region,  Lincoln  Laboratory  sent  an 
analyst  to  Tokyo,  Japan,  during  the  summer  of  1966.  Data  was  collected  from  the 
seismic  tripartite  operated  by  the  University  of  Tokyo's  Earthquake  Research  Institute. 
It  was  hoped  that  this  data  would  provide  a  reference  population  of  small  events  with 
precise  locations  and  magnitudes  so  that  detection  thresholds  and  location  accuracy  of 
LASA  could  be  studied  in  detail,  at  least  for  one  source  region. 

Unfortunately,  the  problems  associated  with  interpretation  of  data  from  the 
tripartite  for  local  events,  particularly  the  magnitude  determination,  were  more  serious 
than  expected.  It  was  possible  to  use  the  data  only  to  determine  the  existence,  approxi¬ 
mate  location  and  origin  time  of  local  events.  The  experiment  was  further  hindered  by 
the  fact  that  there  were  so  many  small  local  events  detected  in  Japan  that  it  was  imprac¬ 
tical  to  save  LASA  tapes  for  each  one.  It  took  several  weeks  to  work  out  acceptable 
communications  and  procedures  for  saving  the  appropriate  LASA  data.  In  the  final 
analysis  there  were  approximately  150  events  that  were  recorded  on  the  tripartite  for 
which  LASA  data  were  available  for  analysis.  This  population  reduced  to  28  events  that 
were  large  enough  to  be  conceivably  detected  by  LASA,  the  remaining  events  being  elim¬ 
inated  when  it  became  evident  that  none  of  the  smaller  events  would  be  detected. 
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For  those  events  that  were  detected  in  off-line  beams,  formed  by  combining 
the  straight  sum  output  of  each  subarray,  the  amplitude  and  period  of  the  best  beam 
for  each  event  was  used  to  compute  the  LASA  magnitude.  The  resulting  magnitudes 
are  not  inconsistent  with  a  LASA  beam  detection  threshold  for  teleseismic  P-waves 
of  about  3.  5. 

The  LASA  locations  were  determined  by  forming  a  grid  of  beams  covering  the 
epicenter  determined  from  local  tripartite  data  and  choosing  the  best  beam.  The  lo¬ 
cation  error  is  magnitude  dependent  at  small  magnitudes.  For  most  events,  the  lo¬ 
cation  indicated  by  the  best  beam  was  within  150  kilometers  of  the  location  determined 
from  local  data,  the  largest  difference  between  the  local  and  LASA  epicenter  being 
400  kilometers. 


II.  DATA  REDUCTION 


During  the  two  months  of  August  and  September  of  1966,  arrival  times  and 
amplitudes  of  P  and  S  waves  were  read  from  recordings  obtained  from  the  three  station 
network  operated  by  the  Earthquake  Research  Institute  in  Tokyo,  Japan.  The  data  ob¬ 
tained  from  the  tripartite  were  similar  in  format  to  data  sent  to  the  U.  S.  Coast  and 
Geodetic  Survey.  The  objective  was  to  reduce  this  data  to  obtain  accurate  epicenters 
and  magnitudes  for  comparison  with  LASA  beamforming  data. 
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Use  of  the  tripartite  to  determine  magnitudes  proved  to  be  a  problem.  The 
time  scale  on  the  chart  recordings  was  too  compressed  to  allow  a  measurement  of  the 
dominant  period,  and  since  the  period  changes  significantly  between  events,  particularly 
at  different  ranges,  it  was  impossible  to  determine  a  consistent  magnification.  The 
instrument  response  for  the  three  stations  in  the  tripartite  is  shown  in  Fig.  1.  The 
magnitude  problem  was  further  complicated  by  the  fact  that  the  instruments  saturated 
at  about  1.  5  microns  so  that  amplitudes  for  large  events  could  not  be  measured. 

The  locations  of  the  three  stations  are  as  follows: 


Tsukuba 

Lat.  : 

36.211°  N 

Long. : 

140. 110°  E 

Dodaira 

Lat. : 

36. 003°  N 

Long. : 

139. 193°  E 

Kiyosumi 

Lat. : 

35. 200°  N 

Long. : 

140.  141°  E 

Data  for  each  event  consisted  of  at  least  two  P-arrival  times  and  at  least  one 
S-P  interval.  Many  events  were  recorded  by  all  three  tripartite  stations.  For  some 
events,  one  station  was  noisier  than  the  others,  and  frequently  provided  only  a  P- 
arrival  time  of  relatively  low  accuracy.  Instead  of  using  a  least -squares  procedure 
for  determining  epicenters,  based  on  a  variable  number  of  measurements,  an  algorithm 
was  used  for  finding  a  unique  pair  of  epicenters  from  two  P-arrivals  and  the  S-P  interval 
from  one  of  these  stations.  The  difference  in  P-arrivals  puts  the  epicenter  (assumed 
to  be  at  the  surface)  on  a  hyperbola,  and  the  S-P  interval  puts  it  on  a  circle,  which 
intersects  the  hyperbola  in  two  points.  The  ambiguity  is  resolved  by  the  third  P- 
arrival,  when  available.  Details  are  given  in  the  Appendix. 
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All  possible  combinations  of  two  P-arrivals  and  one  S-P  interval  were  worked 
up  for  each  event,  and  all  the  resulting  epicenters  for  each  event  were  plotted  on  a 
map.  For  events  having  redundant  data,  the  resulting  epicenters  were  "averaged"  by 
eye,  leaving  one  or  at  most  two  epicenters  per  event.  This  manual  averaging  (which 
is  analogous  to  a  least -squares  procedure)  permitted  the  exercise  of  judgment  in  re¬ 
jecting  solutions  thought  to  be  out  of  line  because  of  poor  data. 

For  each  event,  the  expected  LASA  arrival  time  was  computed  for  one  or  both 
epicenters,  and  the  LASA  data  searched  for  the  event.  Although  LASA  was  not  re¬ 
cording  on  a  24-hour  basis  and  it  was  not  possible  to  save  recordings  of  each  event 
reported,  digital  tapes  were  saved  which  included  the  expected  arrival  times  of  about 
150  of  the  events  located  from  tripartite  data. 

After  the  epicenters  were  computed  for  the  population  of  events,  a  grid  of 
beams  at  50  km  spacings  was  formed  around  the  epicenters.  The  beams  were  formed 
by  phasing  and  adding  the  subarray  straight  sums  using  the  standard  Jeffreys -Bullen 
travel  time  tables  and  station  corrections  prepared  by  Earth  Sciences,  a  Teledyne 
Company.  ^  Beamforming  started  with  the  largest  events  and  worked  down  until  it 
was  evident  that  no  smaller  events  would  be  detected.  In  this  manner,  the  original 
population  was  reduced  to  28  events  that  were  clearly  seen  or  were  considered  to  be 
marginal.  Additional  information  on  these  events  is  given  in  Table  I. 
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TABLE  1 


EVENT  LIST 


1966 

Date 

Origin 

GMT 

Latitude 

Longitude 

Expected 

LAO 

Arrival 

Magnitude 

PDE 

Detected  Card 

8/19 

06  22  19 

32.  1  N 

139.3  E 

06  34  37 

4.8 

Yes 

8/19 

06  26  38 

32.0  N 

139.0  E 

06  38  57 

4.5 

Yes 

8/20 

08  31  45 

19.  1  N 

144.3  E 

08  44  36 

Questionable 

8/20 

09  32  32 

43.  1  N 

140.  6  E 

09  43  44 

5.8 

Yes  61 

9/02 

10  19  55 

37.  0  N 

139.5  E 

10  31  51 

3.  6 

Yes 

9/03 

08  11  39 

43.2  N 

146.5  E 

08  22  44 

4.  6 

Yes  64 

9/10 

02  27  48 

46.6  N 

144. 1  E 

02  38  16 

5.2 

Yes  65 

9/11 

06  51  45 

41.4  N 

143.  0  E 

07  03  13 

4.0 

Yes 

9/11 

17  33  30 

38.  2  N 

137. 3  E 

17  45  28 

3.5 

Questionable 

9/15 

03  27  07 

36.6  N 

138.2  E 

03  39  09 

Questionable 

9/17 

06  32  21 

27.0  N 

142.  5  E 

06  44  49 

3.  6 

Yes 

9/18 

05  22  31 

42.3  N 

142.  8  E 

05  33  51 

5.  1 

Yes  67 

9/18 

08  14  08 

35.  1  N 

142.  7  E 

08  26  02 

3.6 

Questionable 

9/18 

10  28  29 

35.9  N 

139.  8  E 

10  40  29 

3.4 

Questionable 

9/21 

04  02  14 

36.  9  N 

133.  9  E 

04  14  27 

3.6 

Yes 

9/21 

06  03  18 

39.5  N 

144. 3  E 

06  14  48 

3.  6 

Yes 
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9/22 

04  15  31 

37.3  N 

138.6  E 

04  27  26 

4.9 

Yes 

8/20 

09  32  37 

36.  6  N 

138.2  E 

09  44  39 

No 

9/15 

06  19  01 

37. 1  N 

139.5  E 

06  30  56 

No 

9/16 

08  00  01 

37.  1  N 

138. 7  E 

08  11  59 

No 

9/16 

09  49  39 

35.  6  N 

139.  8  E 

10  01  41 

No 

9/16 

10  16  04 

31.9  N 

139. 5  E 

10  28  22 

No 

9/16 

11  36  02 

41.2  N 

143.  4  E 

11  47  29 

No 

9/17 

11  23  37 

39.  1  N 

136. 8  E 

11  35  32 

No 

9/18 

06  42  01 

36.5  N 

140. 3  E 

06  53  57 

No 

9/18 

08  40  53 

39.3  N 

141.  1  E 

08  52  35 

No 

9/21 

10  12  01 

35.  9  N 

142. 7  E 

10  23  53 

No 

9/22 

02  58  41 

32.  1  N 

138. 8  E 

03  11  00 

No 
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III. 


DETECTION  THRESHOLD 


In  an  earlier  report,  it  had  been  estimated  that  the  Montana  LASA,  using 
subarray  outputs  could  detect  50%  of  the  teleseismic  P-wave  arrivals  of  magnitude 
4.  1  or  greater  and  that  therefore  the  corresponding  magnitude  threshold  when  beams 
were  used  should  be  3.  5.  One  of  the  objectives  of  this  experiment  with  Japan  data 
was  to  check  this  assertion.  Figure  2  shows  a  plot  of  the  number  of  events  detected 
versus  magnitude.  The  number  of  events  detectable  on  LASA  beams  is  much  too  small 
to  allow  an  accurate  estimate  of  the  detection  threshold,  but  the  magnitude  of  these 
events  seen  and  the  noise  level  on  the  beams  are  not  inconsistent  with  a  detection 
threshold  of  about  3.  5  for  beams  steered  to  Japan. 


IV.  EPICENTERS  DETERMINED  BY  BEAMSPLITTING 


The  epicenters  obtained  from  the  Japan  tripartite  data  were  thought  to  be  within 
a  LASA  beamwidth  (8°)  of  the  true  epicenters.  It  was  desired  to  test  the  location  ability 
of  the  LASA  by  the  technique  of  "beamsplitting'',  that  is,  careful  determination  of  the 
direction  of  aim  of  the  LASA  beam  that  produces  the  greatest  signal  output.  The  popu¬ 
lation  of  28  events  was  analyzed  in  this  manner.  The  results  produced  12  events  which 
could  be  clearly  seen  and  identified,  five  events  which  were  marginal  and  considered 
questionable,  the  remaining  11  events  being  definitely  not  visible  even  on  filtered  beams. 
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The  locations  of  the  computed  epicenters  were  compared  with  the  locations 


of  the  best  beam  resulting  from  the  grid  of  beams.  The  difference  in  the  two  epicenters 
was  plotted  against  the  LASA  magnitude.  In  most  cases  the  best  beam  was  within  1.  5 
degrees  (166  km)  of  the  computed  epicenters.  The  results  for  the  12  clearly  seen 
events  are  shown  in  Fig.  3.  The  curve  is  a  theoretical  formula  for  beamsplitting 
accuracy  based  on  idealized  models  of  the  signal,  noise  and  the  measured  process. 

The  theoretical  model  of  the  signal  assumes  perfect  knowledge  of  time  delay  station 
corrections. 

By  expanding  the  grid  of  beams  from  the  original  16  to  400  and  measuring  the 

output  power  of  the  signal  in  each  beam  and  then  contouring  the  signal  outputs,  one 

3 

obtains  a  true  picture  of  the  actual  beam  pattern.  Figures  4  and  5  show  the  results 
of  performing  this  operation  on  two  of  the  12  detected  events  listed  in  Table  1. 


V.  CONCLUSION 


Although  the  attempt  to  measure  the  Montana  LASA  teleseismic  detection 
threshold  did  not  prove  definitive,  the  results  are  encouraging.  The  experiment 
should  be  reported  for  other  situations  where  both  the  LASA  and  a  local  network  monitor 
low-magnitude  seismic  activity  in  a  region  distant  from  Montana,  for  example  the 
Kurile  Islands  ocean  bottom  experiment  of  late  1966. 
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The  present  experiment  has  provided  a  meaningful  determination  of  epicenter 
location  error  using  beamsplitting  of  a  single  LASA  for  a  source  region  in  which  the 
station  corrections  are  well  known. 
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APPENDIX 


The  population  of  events  used  in  this  experiment  ranged  in  location  from 
sources  within  the  tripartite  to  epicenters  several  hundred  kilometers  away.  We  have 
no  real  idea  what  the  distribution  of  events  was  in  depth  of  hypocenter,  and  all  events 
were  assumed  to  have  surface  foci.  It  is  therefore  not  to  be  expected  that  a  simple 
algorithm  will  provide  high  accuracy  in  location  and  it  is  also  not  easy  to  infer  from 
known  properties  of  the  region  what  method  to  use  for  epicenter  location. 

The  procedure  we  used  requires  two  basic  formulas:  one  to  give  difference 
in  distance  to  two  stations  from  a  measured  difference  in  P-arrival  times  and  a  second 
to  give  distance  itself  from  a  single  station  using  a  measured  S-P  interval.  If  the 
travel -time  for  P  is  a  linear  function  of  distance,  such  as 

Tp  =  Tq+  (A/Vp)  , 

then  the  first  required  formula  is 

A  -A  =  V  (T  -  T  ) 

2  1  Pv  2  V 

The  measured  P-arrival  times  are  T.  and  T_  and  V  is  assumed  P-speed'.  Note  that 

i.  ^  I 

the  constant,  T^,  is  not  used.  We  used  the  value  7.  5  km/sec  for  V^,  which  is  typical 

4 

of  P^  speeds  in  Japan. 

In  order  to  get  distance  from  the  S-P  interval,  T_  ,  we  used  the  relation 

•bx 


10 


where 


=  ATSP+BTSP 


and 


A  =  7.  86  km/sec 
2 

B  =  0.04  km /sec 


The  small  quadratic  term  was  used  to  increase  the  (A/Tgp)-ratio  from  a  value  near 
8  km/sec,  thought  to  be  typical  of  near  events,  to  a  value  near  10  km/sec  for  a  dis¬ 
tance  of  500  km.  The  value  10  km/sec  is  consistent  with  Jeffreys -Bullen  travel  times 
for  distances  in  the  500  —  1000  km  range.  In  other  words,  we  tried  to  interpolate  be¬ 
tween  a  value  used  in  Japan  for  locating  local  events  and  the  standard  tables  for  events 
at  relatively  great  distance. 

A  limited  amount  of  experimenting  with  other  constants  in  these  formulas 
showed  that  our  events  were  just  as  consistently  located  (i.e. ,  epicenters  from  redundant 
data  on  one  event  were  in  good  agreement)  by  a  linear  relationship  between  distance  and 
S-P  interval.  If  one  assumes  that  both  the  P  and  S  travel  time  curves  are  linear,  and 
that  one  has  S-P  measurements  for  two  stations,  (TQ  )  and  (T  )  ,  as  well  as  their 
P-arrival  times,  P  and  P  ,  then  it  follows  that 


^TSP^2  (Vl 

P2-pi 


a  -  1  , 


where  a  =  Vp/V^  is  the  ratio  of  P  to  S  speeds.  We  measured  this  ratio  for  all  events 
in  the  population,  with  the  average  result  a  =  1.  81.  This  ratio  requires  an  S  speed  of 


11 


4.  14  km/sec,  to  be  consistent  with  our  assumption  of  the  value  7.  5  km/sec  for  P. 

These  speeds  in  turn  imply  the  relation 

A  =  9.2  Tsp 

for  distance  in  terms  of  S-P  interval.  Our  quadrant  formula  predicts  the  ratio  (A/Tgp) 
=  9.  2  at  a  distance  of  300  km.  Since  the  average  distance  of  the  events  from  the  center 
of  the  tripartite  is  several  hundred  kilometers,  our  formulas  for  epicenter  accuracy 
are  in  reasonable  agreement  with  the  data  itself. 
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MAGNITUDE 


Fig.  1.  Instrument  magnification. 
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Fig.  2.  Detected  events. 
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Fig.  3.  Epicenter  error. 
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Fig.  4.  Actual  beam  pattern. 
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Fig.  5.  Actual  beam  pattern. 
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